Alberta Thy 04-10, TUM-EFT 8/10, UB-ECM-PF 10/15, ICCUB-10-029 



Precision determination of r Aj^g from the QCD static energy 



(N 

o 



Nora Brambilla, 1 Xavier Garcia i Torino, 2 Joan Soto, 3 and Antonio Vairo 1 

1 Physik Department, Technische Universitat Munchen, D-85748 Garching, Germany 
2 Department of Physics, University of Alberta, Edmonton, Alberta, Canada T6G 2G^\ 
3 Departament d'Estructura i Constituents de la Materia and Institut de Ciencies del Cosmos, 
Universitat de Barcelona, Diagonal 647, E-08028 Barcelona, Catalonia, Spain 

(Dated: May 29, 2012) 

We use the recently obtained theoretical expression for the complete QCD static energy at next- 
to-next-to-next-to leading-logarithmic accuracy to determine roA^jg- by comparison with available 
lattice data, where ro is the lattice scale and Ajjg- is the QCD scale. We obtain roA^ = 0.637^q'q3q 
for the zero-flavor case. The procedure we describe can be directly used to obtain roA^ in the 
unquenched case, when unquenched lattice data for the static energy at short distances becomes 
available. Using the value of the strong coupling a s as an input, the unquenched result would 
provide a determination of the lattice scale ro. 

PACS numbers: 12.38.Aw, 12.38.Bx, 12.38.Cy, 12.38.Gc 
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The energy between a static quark and a static an- 
tiquark is a fundamental object to understand the be- 
havior of quantum chromo dynamics (QCD) [l[ . Its long- 
distance part encodes the confining dynamics of the the- 
ory while the short-distance part can be calculated to 
high accuracy using perturbative techniques. Perturba- 
tive computations of the short-distance part have been 
performed for many years @, Q and the two-loop cor- 
rections have been known for quite some time now 0- 
6]. When using perturbation theory to calculate the 
short-distance part, the virtual emission of gluons that 
can change the color state of the quark-antiquark pair 
(so-called ultrasoft gluons) produce infrared divergences, 
which induce logarithmic terms, lna s (l/r), in the static 
energy. Those effects, which first appear at the three- 
loop order, were identified in Ref. pj and calculated 
in Ref. [1, Q using an effective field theory framework 
lol 11 1 . That framework also allows for resummation of 
the ultrasoft logarithms [l2| , which may be large at small 
distances r. Very recently, the complete three-loop cor- 
rections to the static energy have become available 
15[ . Combining the results of those calculations with the 
resummation of the ultrasoft logarithms at sub-leading 
order 16|, [l7| , the static energy at next-to-next-to-next- 
to leading-logarithmic (N 3 LL) accuracy, i.e. including 
terms up to order a^ +n ln n a s with n > 0, is now com- 
pletely known. 

In the first part of the letter, we compare the static en- 
ergy at N 3 LL accuracy with lattice data. The comparison 
shows that, after subtracting the leading renormalon sin- 
gularity, perturbation theory reproduces very accurately 
the lattice data at short distances, thus confirming at 



key ingredient to relate low energy hadronic physics with 
high energy collider phenomenology. This constitutes the 
main result of our work. 

The static energy Eo(r) at short distances can be writ- 
ten as 



Eo(r) = V S +A S + S- 



us, 



(1) 



where V s and A s are matching coefficients in potential 
Non-Relativistic QCD (pNRQCD) and <5us contains 
the contributions from ultrasoft gluons. V s corresponds 
to the static potential and A s inherits the residual mass 
term from the Heavy Quark Effective Theory Lagrangian. 
In order to obtain a rapidly converging perturbative se- 
ries for the static potential in the short-distance regime, 
it has been argued that it is necessary to implement a 
scheme that cancels the leading renormalon singularity 
[2f]j |. The use of any such scheme introduces an addi- 
tional dimensional scale (which we call p) , upon which all 
the quantities in Eq. ([1]) depend. We will employ the so- 
called RS scheme [21( , in the same way as it was done in 
Ref. [l6j]. The explicit expressions for Eq at N 3 LL accu- 
racy were presented in Ref. 16] and will not be repeated 
here. We refer to that paper for details. The only new in- 
gredient is that the three-loop gluonic contribution to the 
static potential is now known. At three-loop order the 
static potential presents infrared divergences, which can- 
cel in the physical observable Eq after the inclusion of the 
ultrasoft effects. Therefore, it is necessary to consistently 
use the same scheme to factorize the ultrasoft contribu- 
tions for all the terms in Eq. (JlJ. That way one obtains 
the correct three-loop coefficient for the static energy Eq, 
which is independent of the scheme used to factorize the 

ir 



ultrasoft contributions. Refs. 



15] present the result 



an unprecedented precision level the conclusions reached 

18 ]j| In the second part of ^ or ^ ne purely-gluonic three-loop coefficient of the static 



in previous analyses [16 



the letter, the excellent agreement of perturbation the- 
ory with lattice data allows us to obtain a precise deter- 
mination of the quantity roA^jg, where ro is the lattice 



^ (0) (fol- 
ize again 



potential in momentum space, which we call 
lowing the notation of Ref. 15]). We emphas 
that is scheme dependent. The corresponding coeffi- 
scale and A^jg is the QCD scale (in the MS scheme), a cient in the static energy can be obtained by taking the 
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d-dimensional Fourier transform of the momentum-space 
potential (as calculated in Refs. [14|,|l5j]) and adding to 
it the ultrasoft contribution. In the factorization scheme 
used in Refs. lj, Il5f. the ultrasoft contribution is given 
by Eq. (8) of Ref. |l4| (which we confirm). Note that 
the scheme used in Refs. , [lj| is different from the one 
used in Ref. [I(J. If we then subtract Eq. (34) of Ref. [l6| 
(the ultrasoft contribution in the scheme of that paper) 
from the static energy we get the three-loop gluonic con- 
tribution to the static potential in the scheme of Ref. [3] 
(which we will denote as ra ). By doing that we ob- 
tain 



Co 



a 



(0) 

3, Ref. [Tjl 



/4 3 = 222.703. This was the only 
missing ingredient of the static energy at N 3 LL accuracy 
at the time that paper was written. Note that the value 
of Co above is within the range (215,350) predicted in 
Ref. [3] , and considerably lower than the Pade estimate 
co = 313 commonly used in the literature [Hj]. 

We now compare the perturbative results for the static 
energy with the nt = lattice data of Ref. 23| . All re- 
sults are presented in units of r . As it is explained in 
Ref. flij ]. an appropriate quantity to plot for this com- 
parison is 



E (r) - E a (r min ) + ^ att -(r min ), 



(2) 



where 



is the shortest distance at which lattice data 



is available and r E l att - (r min ) = -1.676 [2j. Now that 
the three-loop static potential is known the normalization 
of the u = 1/2 renormalon singularity i? s , which is nec- 
essary to implement the RS scheme, can be determined 
using one order more in the perturbative expansion of 
the potential. We obtain 



R s 



-1.333 + 0.499 - 0.338 - 0.033 = -1.205, (3) 



which is the value we will use. The rest of the scales 
and parameters are set as in Ref. [lf|. We note that, 
in particular, this means that we have p = i.lbr^ 1 and 



r A^jg = 0.602 [24|. The comparison of the static energy 



with lattice data is presented in Fig. [Ha). We recall that 
all the curves coincide with the lattice point at r = r min 
by construction [as one can see from Eq. @]. We note 
that, due to the singlet-octet mixing in the renormaliza- 
tion group equation for A s (at order r 2 in the multipole 
expansion of pNRQCD), the N 3 LL curve (solid black) de- 
pends on a constant, which we call K<i. Power counting 
tells us that \K^\ ~ ^m$, but the constant is otherwise 
arbitrary. We fix it by a fit to the lattice data, which 
delivers r$K2 = —2.06. The bands in Fig. [Ha) illustrate 
the effect of variations in roA^g = 0.602 ± 0.048 [H 
(yellow lighter band) and the effect of adding the term 
±GVag /r, which is representative of the neglected higher 
order corrections (green darker band); the remaining pa- 
rameters are kept at their original values to obtain the 
bands. The figure confirms that perturbation theory (in 
the RS scheme) reproduces very well the lattice data for 
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FIG. 1. (a) Comparison of the singlet static energy with lat- 
tice data. We plot r a (Eo(r) — Eo(r m i n ) + E, 
function of r/ro and the lattice data of Ref. [2 
The dotted blue curve is at tree level, the dot-dashed ma- 
genta curve is at one loop, the dashed brown curve is at two 
loop plus leading ultrasoft logarithmic resummation and the 
solid black curve is at three loop plus next-to-leading ultrasoft 
logarithmic resummation. The yellow (lighter) band is ob- 
tained by varying roA^jg according to roA^jg = 0.602 ± 0.048 
[24| and the green (darker) band is obtained by adding the 
term ±Cfcz^ /r , for the solid black curve, (b) Same but using 
roA^g- = 0.637iS;gjg (see text). 



the static energy at short distances. We also note that the 
band due to the neglected higher order terms is smaller 
than the one due to the uncertainty in roAj^g-. Those 
facts indicate that we should be able to use the lattice 
data to obtain a more precise determination of roA^r, as 
we describe later. It is worth emphasizing that the aim 
of Fig. Q] is to compare with lattice data for the specific 
r range displayed in the figure and to see if the theoreti- 
cal curves follow the lattice data points starting from the 
point at the shortest distance. The scale p was set to the 
fixed value 3.25^*, which is at the center of the range, 
and it was kept the same for all curves. Moreover, the 
static energy was expressed as a series in a s (p), rather 
than a s (X/r), to reduce the uncertainty associated with 
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R s and the specific implementation of the renormalon 
subtraction. 

We will now assume that perturbation theory by itself 
(after canceling the leading renormalon) is indeed enough 
to accurately describe the lattice data for the range of r 
we are considering, i.e. 0.15ro < r < 0.5ro. That is, we 
are assuming that non-perturbative effects are small and 
can be totally neglected 0- With this assumption, we can 
use the lattice data for the static energy to determine 
roAjjg. We initially consider the set of roA^jg- for which: 
(i) the perturbative series for the static energy appears to 
converge, and (ii) the agreement with lattice is improved 
when increasing the perturbative order of the calculation. 
If we implement the condition (ii) by demanding that 
the reduced \ 2 01 the curves decreases when we increase 
the perturbative order of the calculation, we obtain the 
range (0.58,0.8) for r Aj^g. We now proceed to improve 
the precision of that determination. First we recall that 
the use of the RS scheme introduced the scale p in the 
expressions for the static energy. We have used the value 
p = S.Zbr^ 1 because it corresponds to the inverse of the 
central value of the r range we are comparing with lattice 
data. That way we keep In rp terms from becoming large 
[lfjj], but in principle any p around that value is valid. We 
will exploit this freedom to find a set of p values which are 
"optimal" for the determination of the parameter roAjjg 
by following the procedure we describe next: 



1. We vary p by ±25% around p 
p = 2.44 to p = 4.06@ 



3.25r , i.e. from 



2. For each value of p and at each order in the pertur- 
bative expansion of the static energy, we perform 
a fit to the lattice data. The parameters of the fit 
are njA^jg for the curves from tree level to next- 
to-next-to leading logarithmic (N 2 LL) accuracy (1- 
parameter fits), and roA^jg and rgK 2 for the curve 
at N 3 LL accuracy (2-parameter fit). 

3. We select those p values for which the reduced \ 2 
of the fits decreases when increasing the order of 
the perturbative expansion. 

4. Finally, from the set of p values obtained above, 
we select the ones for which the fitted value of K% 
is compatible with the power counting (we require 
\r K 2 \ < 2 p3). 

The above steps provide us with a certain range of p. We 
then consider the set of fitted values of ^oAjyfg at N 3 LL 



accuracy (denoted as Xi below) for that range. In order 
to give more significance to the better fits, we assign a 
weight to each of the sc,-. We choose those weights to be 
given by the inverse of the reduced \ 2 of the fit. We take 
the weighted average of the Xi as our central value for 
the determination of roAjyjg- and obtain 



x :~ W{Xi = 0.637, 



(4) 



where Wi := Wi/(^2j Wj) and ibi is the weight of the 
point Xi. To estimate the error that we should associate 
to that number, we first consider the weighted standard 
deviation of that set of values and assign it as an error 
to the weighted averag^l, we obtain 



1 



i-E 



3 W ) 



^2wi(xi -x) 2 =0.012. (5) 



Fig. [2] shows the obtained fit values of rnAjjg at N 3 LL 
accuracy and at N 2 LL accuracy, for the different values 
of p. The size of each point in the plot is proportional 
to its weight. Additionally, we also consider the dif- 
ference between the weighted averages computed using 
the N 3 LL result and the N 2 LL result and assign it as a 
second error to the weighted average at N 3 LL accuracy 
(we compute the weighted average at N 2 LL accuracy in 
the p range obtained after step 3 of the procedure de- 
scribed above, since step 4 applies only to the N 3 LL re- 
sult). The value of that difference is 0.011. We sum 
those two errors linearly and obtain 0.637 ± 0.023. We 
present in Table U the obtained values for rrjAjgg- follow- 
ing this procedure at different orders of accuracy. We 
emphasize that the error assigned to the result must ac- 
count for the uncertainties associated to the neglected 
higher order terms in the perturbative expansion of the 
static energy; in that sense, we note that Table Q] shows 
that our errors at a certain perturbative order always in- 
clude the central value at the next order, which gives us 
confidence in the reliability of the procedure. In order to 
further assess the systematic errors steming from our pro- 
cedure, we have redone the analysis using two additional 
weight assignements: (1) p- value weights, this analysis 
gives ^oAf^g- = 0.638 ± 0.031, and (2) constant weights, 
this analysis gives njA^g- = 0.634 ± 0.028. These two 
numbers are compatible with our previous value. In our 
final result, we will quote an error that covers the whole 
range spanned by the three analyses, which we consider 
a fairly conservative estimate. Finally, we mention that 
there is also an error associated to each of the fit values 
Xi coming from the variance of the lattice data pointfl 



1 The leading genuine non-perturbative contribution, which is pro- 
portional to the gluon condensate, is of order Ai^r 3 /a s (l/r), 
and hence parametrically suppressed according to the counting 
of Ref. [HI, which assumes Aj^g- ~ Q s /^. 

2 We have used steps of 1.6 X lO -3 ^ 1 to do that. 



3 A similarly motivated procedure to estimate theoretical errors 
has been used, for example, in Ref. [25| , 

4 To obtain this error we consider 1-parameter fits, with Ks fixed. 
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FIG. 2. Fit values of roA^g- for different values of p. The 
range of p displayed in the plot corresponds to that obtained 
after step 3 (N 2 LL points) or step 4 (N 3 LL points) of the 
procedure described in the text. The size of each point is 
proportional to its assigned weight. 



Accuracy 


r oA w 


next-to- leading 


0.946 ± 0.039 ± 0.80 


N 2 LL 


0.648 ± 0.006 ± 0.30 


N 3 LL (no p.c.c.) 


0.619 ±0.04 ±0.029 


N 3 LL 


0.637 ±0.012 ±0.011 



TABLE I. Values of foA^s obtained at different levels of ac- 
curacy. N 3 LL (no p.c.c.) stands for N 3 LL accuracy without 
imposing the power counting constraint on K%, i.e. omitting 
step 4 in the procedure described in the text. The first er- 
ror corresponds to the weighted standard deviation and the 
second one to the difference with the previous order. 



The error induced in this way in x is much smaller than 
the other ones we are considering (due to the fact that 
the lattice points have very small error bars) and can be 
neglected. 

According to the discussion above, our final result for 



+0.032 
0.030- 



MS 



(6) 



0.602 ± 



Our result is compatible with the value roA 
0.048 given in Ref. .241], which we had been using previ- 
ously, but has a smaller error. We also mention that a 
new (preliminary) lattice calculation, which determines 
r o A^Js from the ghost and gluon propagators, obtains a 
result, njA^jg = 0.62 ± 0.01, which is similar to ours [2fj| . 
We present in Fig. QJb) a comparison of the static en- 
ergy with lattice data using the value obtained in Eq. ([6]) 
for roAj^g, and the best fit value for roK 2 when that 
value of ^oAjyjg is used, namely r K 2 = —1.39. We note 
that the error bands due to uncertainties in roA^g and 
higher-order terms in Fig. QJb) are of comparable size. 

We would like to emphasize that exactly the same pro- 
cedure we have described here could be used in the un- 
quenched case. If unquenched lattice data for the static 



energy at short distances were available we could obtain 
an unquenched value for ^ciAms" from it. Combining that 
result with the value of the strong coupling a s , deter- 
mined from other sources, would provide a determina- 
tion of the lattice scale ro . This determination would be 
model independent and alternative to other determina- 
tions like, for instance, the one using the 1S-2S bottomo- 
nium radial excitation energy 27 1. 

In summary, we have updated the comparison of the 
static energy at N 3 LL accuracy with lattice data by in- 
cluding the recently calculated purely-gluonic three-loop 
contribution to the static potential, which was the last 
missing ingredient, (see Fig. [IJ. We confirmed that, af- 
ter canceling the leading renormalon singularity, pertur- 
bation theory can accurately reproduce the lattice data 
at short distances. By taking advantage of this fact, we 
have obtained a new determination of roAirro-. Our result 



for the zero-flavor case is ^oA-p^g = 0.637^q'q3q, which 
improves the precision of the N 2 LL determination (see 
table HJ) by an order of magnitude. 
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